A suggested 4x4 structure in underdoped cuprate superconductors: a Wigner 

supersolid. 
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A wave function is proposed for the "4 x 4" inhomogeneous structures observed on cuprate super- 
conductors. It is based on the Gutzwiller-RVB technique proposed in recent papers, and consists of 
a Wigner sohd of hole pairs embedded in a sea of d-wave spin singlet pairs. Arguments are given 
that the nodal quasiparticles may remain unscattered and even superconducting on such a structure. 



A number of recent STM experiments on underdoped 
cuprate superconductors have shown evidence of struc- 
ture with a Bravais lattice close to 4a x 4a in the Cu02 
planes. A possibly similar structure has been observed 
in vortex cores in BSSCO, with a fractionally larger lat- 
tice constant; and we also call attention to the often- 
noticed anomalies of T^ and other parameters in LSCO 
at a doping of a; = i which to my knowledge have never 
been satisfactorily explained, and may be caused by the 
same structure. Both the tunnelling spectrum and the 
I structure seem often to indicate that there is super- 
conductivity at low temperatures in materials with these 
structures. 

I propose here a microscopic description of these phases 
that is not inconsistent with their arising at low dop- 
ing levels in an RVB supcrconductor-Mott insulator sys- 
tem. I am suggesting that they are in essence a two- 
dimensional crystal of singlet d-wave pairs of holes, com- 
mensurate in many cases with the underlying lattice, and 
existing within a background of a d-wave RVB of singlet 
pairs. (Which latter is our model of the electronic nature 
of the pseudogap state in these materials.) Our model 
differs in essential ways from those of Chen et al. jll and 
Lee et al. [31, but has in common that we all propose a 
crystal of holes. The model of Lee et al. explicitly re- 
jects hole pairs in favor of individual holes, but seems to 
us to fail to explain why the structure appears when the 
material is still superconducting and still exhibits quasi- 
particle nodes-clearly the pairs are a reality, in spite of 
the Coulomb repulsion between the two holes. It also 
loses the feature of explaining the "i" phenomenon and 
of why the doping level at which these observations occur 
is never as low as i . The model of Chen et al. is based on 
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a proposed symmetry between antiferromagnetism and 
d-wave superconductivity which many consider problem- 
atic, and sees the crystal as a "Wigner crystal", which 
to my mind implies an insulating state. I should note 
that in a number of numerical simulations on underdoped 
multi-leg ladders (White and Scalalpino Q), 4x4 squares 
centered on a hole pair in a single plaquette are observed 
as one of a number of inhomogeneous states. 

In some recent papers 0, Q the author and collab- 
orators have returned to the early insight p] that the 
ground state of the Cu02 planes in the cuprate supercon- 
ductors can be modelled as a Gutzwiller projection of a 
d-wave BCS superconducting state. Our point of view is 



to construct an effective Hamiltonian that operates only 
within the manifold of "lower Hubbard band" states - in 
first approximation, the "t-J" Hamiltonian - and then 
to recognize that the eigenstates of such a Hamiltonian 
must be general Gutzwiller projected states. We treat the 
states before projection by a Hartree-Fock-BCS approx- 
imation, that is we find variationally the best product 
of one-electron functions possible. As in Hartree-Fock, 
the mean field equations which arise from the variational 
procedure also specify the quasiparticle excitations and 
their energies - that is, there is effectively a Koopman's 
theorem for this system. 

In studying point-contact tunneling with this tech- 
nique '&\ we found it useful to use a formulation of the 
ground state wave-function which to our knowledge was 
first given by Laughlin Q: 

I*) = exp{iS)pG X [Z]"— x \<^BCs), (1) 

where Upairs is the number of hole pairs. 
Here, the Gutzwiller projector 

Pg ^ lUil - n,,^n,,^] (2) 

i 

is a projector that removes all doubly-occupied sites but 
otherwise leaves amplitude and phase relations unaltered 
in the d-wave BCS wave function \^bcs)- The product 
in Eq. [3 is over all sites. The BCS function is that ap- 
propriate to the Fermi level for 1 — x electrons, where 
X is the doping fraction. |^) is assumed to be deter- 
mined variationally as in Ref. 4]. Finally, the canonical 
transformation exjp(—iS) transforms the true Hamilto- 
nian into the projected form of the t-J Hamiltonian, and 
correspondingly its inverse transforms the projected wave 
function into the true one. 

Laughlin's innovation [3 was to make explicit the "fu- 
gacity factor" Z raised to the number of hole pairs, which 
is necessary to adjust the populations of the electron 
states on a given site. It is easily calculated that in a 
uniform state Z = 2a; /(I -I- a;), with x the doping frac- 
tion. In Refs. [3, 0, lalthis population adjustment is 
done by fiat. In Ref. [^ the Z factors were, further, 
incorporated into the BCS function, which modifies the 
definition of quasiparticle excitations in such a way as to 
give good agreement with tunneling measurements, but 
we emphasize that the wave function assumed there is 



just Eq. n This variant of the wave function makes 
clear that Z plays the role of a Bose condensate ampli- 
tude - though we should emphasize that our theory does 
not involve assuming a "holon condensate." 



Here, 7(k) = cosik^ — k'^) + cos{ky — k' ) and 
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FIG. 1: Proposed 4x4 structure. Plaquets with crosses are 
nominal sites of hole pairs. The structure may also be thought 
of as "columnar" valence bond structure commensurate with 
hole-pair "liquid crystal". 

What is being proposed here is that we approximate 
the wave function for a density wave or Wigner-like crys- 
tal as described above not primarily by changing the BCS 
function in Eq. ^but by requiring the fugacity Z to vary 
from site to site. Its logarithm is the chemical potential 
for holes and thus must track the Madelung potential of 
the superlattice. The average fugacity is arranged to lead 
to the correct doping level overall, but (for instance) Z 
may be taken to have one value, Zi, on the central four 
sites of a 4x4 square plaquette, and another, smaller 
value Z2 on the remaining 12 sites (or, if desired, one 
may specify still a third value for the corners.) Then the 
hole density will be a maximum on the central plaquette, 
and will be small in the intermediate regions. The con- 
trast in Z values cannot in the end be very great, because 
Zmax will likely not be greater than the value at optimal 
Tc , about i (see Fig. 1). 

The gap equations Eqs. Q and |21 which result from 
variation of |^) are, in the renormalizcd mean- field the- 
ory approximation of Ref. Q , 



A{k) = Ji2-ZfY. 



A(kO 



e(k)=Z6(k) + J(2-Z)2^7(k) 



e(ko 

E{k') 



(3) 



(4) 



Eik) = V?(k)2 + A(k)2 



(5) 



is the standard BCS expression. The important thing to 
note is that the true kinetic energy e(k) is renormalizcd 
relatively by a factor Z/{2 — Z)^ relative to J. The under- 
doped regime is defined by Zt ^ 4J, so that over most 
of the region near the Fermi surface, i?(k) is relatively 
weakly dependent on Z . Thus to a zeroth-order approx- 
imation the gap equations Eq. |21are not affected by pe- 
riodic spatial variation of Z, justifying our basic ansatz 
that the function before projection may be assumed not 
much changed. This signals the important fact that the 
CDW lives in an RVB background, not in a conventional 
band or in an antiferromagnetically (or even spin-glass) 

ordered state. The fact that ^ of a hole relative to the 
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Mott-Hubbard insulator corresponds to a unit cell of 16 
sites does not add up in conventional band theory. 

Of course, the periodic variation of Z will cause Bragg 
scattering of the quasiparticles, which can be repre- 
sented by dividing the spectrum up into 16 sub-bands 
and describing the perturbation as a matrix in band 
indices which opens Bragg-scattering gaps at the sub- 
zone boundaries. I don't feel that Fermi surface nesting 
plays much of a role here. Treating the resulting self- 
consistency problem accurately is beyond the patience 
or ability of the present author to solve directly. But 
there are a number of insights one can recognize. First, 
there will only be one of these subbands which contains 
the nodal quasiparticles, and the lowest part of the spec- 
trum around the nodes will be relatively little scattered, 
since the density of final states goes to zero as E'^ at low 
E. The experimental observation that the nodes sur- 
vive and are coherent in the density-wave state is thus 
confirmed. Second, we can expect the spectrum near 
the antinodes to be severely broken up into rather flat 
"optical" bands, because both the gap and the under- 
lying kinetic energy are rather flat in this neighborhood. 
The superfluid stiffness coming from these portions of the 
zone will be severely reduced by this scattering, and these 
excitations may move diffusively rather than coherently. 

Either from this insight, or on general principles, we 
can expect that the superfluid stiffness ps, which in the 
pure case is known to be proportional to Z, will be re- 
duced more by the low values of Z between the peaks 
than by the higher value at the peaks, and since Tc is pro- 
portional to superfluid stiffness [j], Tc will be degraded. 
The nodal quasiparticles, which are responsible for the 
temperature dependence of p^, are little changed. 

It remains to discuss the motivation for this structure. 
It is obviously the long-range Coulomb interaction that 
furnishes the energy gain, and the stiffness of the hole 
wave function which opposes the deformation. We can 
suppose that Z is unlikely to exceed its value for opti- 
mal doping of 20%, about 0.33, on the central plaquette, 
so that the decrease on the periphery is not severe - I 
estimate 0.2 or so, or a; = 0.1. 



The energy gain is something like the Madelung en- 
ergy of the charge distribution of the "hquid crystal" 
which will be of the order of the square of the charge 
contrast divided by the lattice constant (and corrected 
by the dielectric constant). Thus it scales as 



Ecoulomb OC [SZ)"^ X 1/d 



(6) 



where d is the super lattice constant. 

The excess kinetic energy, on the other hand, will con- 
tain a gradient squared which is proportional to {5Z/dY . 
Clearly there is some d large enough so that the Coulomb 
energy gain wins. It is necessary that d^ = 2/x. Other- 
wise, the Madelung estimate is not correct, and we would 
have to include the Coulomb repulsion of the extra pairs. 
Numerical estimates show that it is reasonable that the 
opposite variations of Coulomb and kinetic energy would 
match for d = Aa. 

The above is not a complete theory. In particular, 
we have not calculated the energetics explicitly. But it 
is based on an explicit wave function, so that numeri- 
cal calculations can be carried out when desired. It does 
have the advantage that it leads to an explicitly supercon- 
ducting state closely related to the uniform state which 
explains many of the properties of the superconductors. 
There are several further questions which remain. 

First, we need a reason why the lattice constant seems 
often to be slightly larger than 4.0, according to Fourier 



transforms of the tunnelling current. The only real-space 
images of the microscopic structure (Davis [2|) show how- 
ever that the structure is made up of domains which are 
not continuously connected. The packing of such a gran- 
ular structure will not be perfect and will have a smaller 
average wave number. We can speculate that commen- 
surability with the lattice favors a definite hole density - 
^, we suppose - and that the grain boudaries adjust the 
net charge. Thus the deviation from 4a provides sugges- 
tive evidence for, rather than against, the lattice of hole 
pairs. 

A second puzzle can be resolved in much the same 
way. The vortex core state at optimal doping exhibits 
a similar modulation, albeit with a lattice constant even 
larger than the previous one - about 4.5a is quoted. In 
the vortex core the superfluid stiffness is reduced by the 
large supercurrent. I expect the amplitude stiffness to 
gradients of Z and the phase stiffness which causes su- 
percurrents to be the same, and if one is reduced the 
other will be, and again the Coulomb energy will be fa- 
vored relatively to the kinetic energy. But here we expect 
the doping to be farther from g and the grains to be even 
smaller, hence the larger deviation from 4.0. 

I would like to acknowledge the sharing of his data 
and his interpretations thereof by J. C. Davis, and useful 
discussions with N. P. Ong. 
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